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ABSTRACT 

Recent 21 cm radio observations of H/ regions in the Small Magellanic Cloud, 
have revealed spatial power spectra of the intensity, which are quite similar in shape 
to those previously deduced for the Galaxy. The similarity, in spite the differences in 
the physical parameters between the Galaxy and the SMC, suggests that the shape 
of the power spectra reflects some underlying mechanism which is not too sensitive 
to the environmental specifics. In this paper we present an interpretation for the 
observational power spectra in terms of a large scale turbulence in the interstellar 
medium, in which the emitting H/ regions are embedded. The turbulence gives rise 
to density fluctuations which lead to the observed intensity fluctuations, in the H/ 
regions. The observational power spectra are used to deduce the turbulence spectral 
function. In the SMC, the turbulence largest eddies are comparable in scale to the 
SMC itself. This implies that turbulent mixing should have smoothed out any large 
scale abundance gradients. Indeed, this seems to be the case, observationally. 

The turbulence is also expected to amplify and shape up the large scale magnetic 
field. Indeed, the observational data indicate the existence of a large scale disordered 
field of the strength expected from energy equilibrium with the turbulent velocity field. 
The large scale turbulence is most probably generated by instabilities in the large 
scale flows induced by the tidal close encounter with the LMC ~ 2 x 10 8 yr ago. The 
life-time of the largest eddies is ~ 4 x 10 8 yr so the turbulence had not yet enough time 
to decay and persists even though the energy source is no longer there. 

Subject headings: ISM: structure - turbulence - ISM: kinematics and 
dynamics - ISM: abundances - ISM: magnetic field 

1. Introduction 

Crovisier & Dickey (1983), and Green (1993) have used 21 cm radio observations 
to derive the intensity spatial power spectra for our Galaxy. Recently, Stanimirovic et 
al. (1999) obtained Hj spatial power spectra for the Small Magellanic Cloud (SMC). 
Some of the observations were interferometric and yielded the power spectrum directly, 
and some single dish observations from which the power spectrum was calculated by 
Fourier transform. 
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The resulting power spectra exhibit a power-law behavior, p(q) oc g 7 , over a wide 
range of spatial scales: ~ 10-200pc in the Galaxy, and ~ 30pc-4kpc for the SMC. Here 
q is the absolute value of the 2-dimensional wavenumber in the plane of the sky, and 
the observational power-law indices are 7 ~ —(2.8 -7- 3) for the Galaxy and 7 ~ — 3 for 



Given the differences in the spatial scales for which the spectra were obtained, 
and the differences in the physical parameters between the Galaxy and the SMC, the 
similarity of the indices is remarkable. It suggests the existence of some underlying 
mechanism that is not sensitive to the specific physical parameters of the environment. 
Stanimirovic el al. (1999) attributed the power-law shape to a fractal structure of H/ 
clouds hierarchy. 

The observational power spectra reveal the existence of intensity correlations over 
a wide range of spatial scales, and this naturally suggests that turbulence in the 
interstellar medium (ISM) is the sought for underlying mechanism. Such an approach 
was suggested by Lazarian (1995) and followed up by Lazarian & Pogosyan (2000). In 
this paper we present an interpretation of the observational power spectra in terms 
of density fluctuations which are generated by a large scale turbulence in the ISM, in 
which the emitting neutral hydrogen regions are embedded. The model determines 
the power spectrum for any given spectral function of the underlying large scale 
turbulence. Therefore, the observational power spectra can be used to study the above 
large scale turbulence. From the observational power spectra, we deduce the spectral 
function of the underlying turbulence: F(k) oc k~ 2 , with k denoting the absolute 
value of the 3-dimensional wavenumber. As detailed in the Discussion section, this 
spectral function is quite common in various components of the ISM and characterizes 
turbulence in a compressible medium. 



The observed lines of sight are generally optically thin in the 21 cm line. For the 
SMC, Stanimirovic et al. (1999) note that self absorption effects are limited to narrow 
ranges in velocities and directions. They find that for most of the covered area the 
emission is from regions with column densities, integrated over the entire velocity span, 
of Nfjj < 2.5 x 10 21 cm -2 . Each of the power spectra corresponds to a velocity range 
which equals 1/6 of the entire velocity span. Thus, one can estimate that typically 
each power spectrum corresponds to a column density Njjj ~ 4 x 10 20 cm~ 2 . 

For the outer arms of the Galaxy, for which Green (1993) obtained the best fitting 
power spectra, we can estimate the column density using the relation (Spitzer 1978) 



where Tb and v are in K and km/s, respectively. For the outer arms, the average 



the SMC. 



2. 



The Model Outline 
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Tb ~ 15K over velocity groups of 3 channels, each of width 2.64km/s (Green 1993), 
yields N Hl ~ 2 x 10 20 cm" 2 . 

To estimate the physical depth D of the emitting neutral hydrogen region, 
contributing to each power spectrum, we take a number density of neutral hydrogen 
of n = 10cm -3 . This value is intermediate between n ~ 5cm~ 3 corresponding to warm 
neutral hydrogen with temperature T <; 200K , and n ^ 30cm -3 of the cool neutral 
hydrogen with T <J 100K (Wolfire et al. 1995). This translates for the Galaxy outer 
spiral arms to a physical depth D ~ 7pc and for the SMC to D ~ 15pc. 

The power spectra were grouped by the observers in radial velocity intervals. 
It is interesting that, both for the Galaxy and the SMC, the power spectra for the 
different velocity groups are similar. In the Galaxy, adopting a rotation curve to 
deduce distances, each such velocity group represents a range of distances. Thus, in 
Green (1993) this range for the outermost spiral arms is ~ 2.2kpc while for the Perseus 
data the range is ~ 0.6kpc. Each such distance range is much larger than the above 
estimate of D ~ 7pc. Therefore, it follows that the intensity is contributed by disjoint 
subregions of neutral hydrogen with total depth equaling D, with each subregion 
having a depth which is only a fraction of D, e. g. ~ lpc. 

In the SMC, the total velocity span probably represents large scale flows, of few 
kpc, that are the result of tidal interactions with the LMC and with the Galaxy. 
Indeed, N-body simulations by Gardiner & Noguch (1996) indicate a very strong tidal 
impact imparted to the SMC by the LMC about ~ 2 x 10 8 yr ago, with a duration of 
~ 1 x 10 8 yr. The resulting simulated velocity field exhibits velocities in the range of 
tens of km/s. This is consistent with the observational velocity span of ~ 80km/s of 
the neutral hydrogen emission as reported by Stanimirovic et al. (1999). The velocity 
range of each of the power spectra, corresponds to a distance interval ~ lkpc which is 
much larger than the total depth of the emitting regions, D ~ I5pc contributing to 
each power spectrum. Therefore in this case too, the emission arises from a number of 
disjoint zones, each of depth which is a fraction of D, e. g. ~ 2pc. 

We note that in both cases, the depth of each of the disjoint emitting regions is 
also much smaller than the scales, in the plane of the sky, probed in the power spectra: 
for the Galaxy up to ~ 200 pc and for the SMC up to 4 kpc. Therefore, the emitting 
H/ regions occupy a small fraction of the ISM volume, in which they are embedded. A 
small space filling factor of neutral hydrogen is implied also in the fractal interpretation 
of Stanimirovic et al. (1999). 

In the optically thin case, the intensity is proportional to the column density of 
neutral hydrogen in the line of sight (Spitzer 1978). Thus, the intensity fluctuations 
are due to density fluctuations in the disjoint neutral hydrogen regions. In the 
present model, the density fluctuations are regarded to be a response to a large 
scale 3-dimensional turbulence in the ISM, in which the neutral hydrogen regions are 
embedded. 
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Since the observational optical depths are small, the following picture suggests itself 
: the Hi subregions (clouds) are waved by a large scale turbulence in the embedding 
medium. In this picture the neutral hydrogen regions are regarded as "passive 
markers" which follow the turbulence field but do not feed back on it dynamically. 
The turbulence induces density fluctuations in the neutral hydrogen subregions that 
manifest as intensity fluctuations. 

To be specific, we consider a simple model of density fluctuations in thin slabs 
or sheets parallel to the plane of the sky, which are separated along the line of sight 
direction by distances much larger than their depth. This geometry is consistent with 
the above estimates of the depth of the emitting subregions. Hi regions are generally 
assumed to have a filamentary structure and filaments are expected to have random 
orientations. This is not to say that all Hi emitting regions are necessarily filamentary. 
However, in our case the contributions to the power spectrum from a given angular 
separation in the plane of the sky will predominantly be those from regions at roughly 
the same distance along the line of sight. So, instead of summing over the actual 
filaments, contributing to each power spectrum, we effectively sum over thin slabs 
parallel to the plane of the sky and separated by distances large compared to their 
depth. One may regard the slab family as mathematically spanning the contributions 
of the actual physical filaments. 

For low galactic latitudes, where most of the neutral hydrogen is located, there 
may indeed be a preference for filaments alignment in the above thin slab geometry, as 
a result of the galactic differential rotation. In the SMC the rotation with respect to 
the Galaxy is expected to have a smaller effect. 

In the next section, we show that the observational power spectra are essentially 
the sum of the spectra contributed by each such thin slab. The mutual contributions 
coming from different slabs are much smaller. 



We wish to relate the observational two-dimensional power spectrum with the 
underlying three-dimensional turbulence. This will provide an interpretation of 
the observational power spectrum and at the same time could hopefully yield clues 
regarding the characteristics of the underlying 3-dimensional turbulence. 

The observational 2-dimensional power spectrum is defined as the Fourier 
transform, over the plane of the sky, of the autocorrelation of the fluctuating brightness 
temperature 6b = Tb— < Tb >• Here Tb is the brightness temperature and < Tg > 
its average over the plane of the sky. Thus, the power spectrum p{q) as function of the 
2-dimensional wavenumber in the plane of the sky, q, is 



3. Hi Spatial Power Spectra and ISM Turbulence 
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Here R, and r are in the plane of the sky, and < 9B(R)9B(R+r) > is the autocorrelation 
of 9b obtained by averaging over all positions R in the plane of the sky. Assuming 
homogeneity and isotropy in this plane, the power spectrum depends only on the 
absolute value of the wavenumber, so that p(q) = p(q) with q = \q\. In the optically 
thin case, the brightness temperature is proportional to the column density of neutral 
hydrogen in the line of sight (Spitzer 1978), thus in the present case 

T B {R) = BN Ul (R) = BJ2 [ dl+Ll/2 n(R,z)dz (3) 

• Jdi-Li/2 

with n the number density and B a constant. The summation is over the disjoint thin 
slabs with depth Lj, located at distances dj. As mentioned before, each depth Lj is 
much smaller than the slab extent on the plan of the sky and also much smaller than 
the distances between slabs. The fluctuating part of the brightness temperature is thus 
expressed as 

9 B (R) = BJ2 [ dt+LJ2 n(R, z)dz (4) 

where n = n— < n > is the fluctuating number density, and < n > is the average 
number density in the slab. The autocorrelation of 9 B is thus 

<9 B {R)9 B {R + r) >= B 2 J2 / <h(R,z')n(R + r,z) > dzdz' (5) 

where, again, the angular brackets represent averaging over all positions R in the plane 
of the sky. 

From Crovisier & Dickey (1983), Green (1993), and Stanimirovic et al. (1999) it can 
be estimated that, typically, the root mean square of 9 B , 9 Brms ~ (0.1 0.2) < Tb > 
implying that the root mean square of n, n rms ~ (0.1 -j- 0.2) < n >. Thus, although 
the medium is compressible, the density fluctuations giving rise to the brightness 
temperature fluctuations are not large. This renders plausible the assumption that the 
density fluctuations are driven by a velocity turbulence and are actually a "passive 
scalar" , i. e. trace the velocity fluctuations but do not feed back upon them. In this 
case, the 2 -point correlation of a density fluctuations field driven by a turbulent 
velocity field, is proportional to the 2-point correlation of the latter (see e.g. Lesieur 
1997) 

< h(R, z')h(R + r, z) >= A 2 < v(R, z') ■ v(R + r,z)> (6) 

with A a constant. For simplicity we further assume that the underlying 3-dimensional 
turbulence is homogeneous and isotropic so that the plane average equals the ensemble 
average 
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< v(R, z') ■ v(R + f, z) >= C 3 (r, z - z') (7) 

where C3 is the 2-point correlation of the underlying 3-elimensional velocity turbulence, 
defined in terms of ensemble average 

C 3 (f) =< v{?) ■ v(f+ r>) >= J ^(ky^ifk (8) 

In equation (8), unlike in the previous equations, the position vectors are 3-dimensional 
and the brackets denote ensemble average. C3 is expressed in terms of the 3-dimensional 
spectrum <3?(fc) which in the homogeneous and isotropic case depends on k, the absolute 
value of the 3-dimensional wavenumber k, 

m = m = . o) 

Here, F(k) is the spectral function of the 3-dimensional velocity turbulence. It equals 
twice the kinetic energy per unit mass and unit wavenumber in the turbulence. 
Accordingly, the mean squared value of the turbulent velocity is 

00 

F(k)dk (10) 

ko 

with irky 1 being the largest scale in the turbulence. 

We can now go back to equation (2) and use equations (5), (6), (8) and (9) to 
express p(q) in terms of F(k) 



= E / / / <%'> k *Y q - r - iq - r+lk ^- z >dfc, d 2 q> d 2 r dz dz' , 

' ' j~j Jdi—Li/2 Jdj-Lj/2 J 

(11) 



leading to 

i+Li/2 rdj+Lj/2 



®(q,k z ) V / / e ik ^ z - z,) dzdz' dk z = (12) 

\i1 ' Jdi-Li/2 Jdj-Lj/2 j 



dk z 



Since \di — dj\ » Li,Lj, the contribution of the second term which results from 
inter-slab correlations is much smaller than that of the first term that is due to 
correlations within each slab. Keeping only the first term gives 
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The assumed homogeneity and isotropy in the plane implies that p(q) = p(q), and use 
of equation (9) leads to 



t A 2 B 2 _ t2 rooF L/q 2 + k 2 ) sa?(k z Li/2) „ 

with q being the absolute value of the 2-dimensional in-plane wavenumber, q. Here 
k z denotes the wavenumber in the z direction, perpendicular to the plane of the sky, 
so that k 2 = q 2 + k 2 . The power spectrum is thus the sum of the power spectra due 
to each of the thin slabs. Equation (14) determines the power spectrum p(q) for any 
given spectral function F(k) of the underlying 3-dimensional turbulence . 



4. Application to the Observational Power Spectra 

It is of interest to use equation (14) in order to deduce information about F(k) 
from the observational power spectrum p{q). In the general case this cannot be done 
uniquely. Nevertheless, for power-law functions it is possible. 

Since the depth of each slab is much smaller than the observed scales in the plane of 
the sky for which the power spectra were obtained, the wavenumbers satisfy q « 7r/Lj. 
From equation (14) follows that in this case, for an underlying 3-dimensional turbulence 
spectrum which is a power- law, F{k) oc k~ rn , the 2-dimensional power spectrum is 
also a power-law, p(q) oc q~ m ~ l . To see this, note that each of the integrals in equation 
(14) is now of the form 



where 



{k z U/2f 
-m/2-i sin 2 (yqLi/2) 



1{q)= l + (yqU/2) 2 dy (16) 

and y = k z /q. 

A plot of I(q) is presented in figure 1 for three representative values of m = 1, 2, 3. 
The wavenumber q is in units of qo, corresponding to the largest scale in the plane of 
the sky. A value qoLi/n = 0.005 typical to the Galaxy was adopted; for the SMC a 
typical value is 0.002. For each m, I(q) is essentially constant down to scales in the 
plane of the sky that are ~ 100 times smaller than the largest scale. Only for very 
small scales, which turn to be smaller than the observed ones, I(q) oc q^ 1 . For the 
SMC value of qoLi/ir, the range of scales for which I(q) is constant, would be wider. 
Thus, each of the terms in equation (14) is a power-law with the same index — (m + 1). 

The above results can be understood qualitatively quite easily. The term 
sm 2 (k z Li/2)/(k z Li/2) 2 effectively limits the integration over k z to the range 
~ 4- n /Li. Approximating it by its average value, yields 
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%)«y o (i + y 2 ) 



-m/2-1 



(17) 



Indeed, in the limit of q << tt/Li which applies here, I(q) is constant. 

The power spectrum for each of the velocity groups is a sum of the spectra of the 
individual slabs. Since each of these power spectra is a power-law with the same index, 
the combined power spectrum will also be a power-law with this index. We conclude 
that a power spectrum p(q) oc g™ m_1 i s the result of an underlying 3-dimensional 
turbulence characterized by a spectral function F(k) oc k~ m . The observational values 
of the index of the power spectrum, imply that m ~ 1.8 -j- 2 for the Galaxy and m ~ 2 
for the SMC. 
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Fig. 1. — Log /(g) as function of Log q, (in units of qo), for qoLi/ir = 0.005. The curves correspond 
to m = 1,2,3 in this order, with m = 1 being the uppermost one. 



5. Discussion 

To interpret the observational spatial power spectra we have assumed that the 
density fluctuations that produce them are a response to a large scale 3-dimensional 
turbulence in the ISM, in which the emitting neutral hydrogen regions are embedded. 
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We derived an expression for the intensity power spectrum in terms of the underlying 
3-dimensional large scale turbulence. 

The power spectra were seen to be the sum of contributions from thin slabs of 
neutral hydrogen, parallel to the plane of the sky, with depth much smaller than their 
extent in the plane of the sky. The depth is also much smaller than the distances, along 
the line of sight, between slabs. In this case, we found that a power-law 3-dimensional 
turbulence spectrum, F(k) oc k~ m , leads to a power-law 2-dimensional intensity power 
spectrum, p(q) oc q~ m ~ 1 . From the observational values of the power spectrum indices, 
it follows that m ~ 2 for the SMC and m = 1.8 — 2 for the Galaxy. 

5.1. Turbulence spectrum 

This 3-dimensional turbulence spectrum resembles the Kolmogorov spectrum 
characterized by m = 5/3 but differs in the value of the index which is closer to m = 2. 
For the Kolmogorov spectrum the turbulent velocity on scale I satisfies Vt(l) oc I 1 / 3 , 
while in the present case vt{l) oc I 1 / 2 . A velocity-size relation with index ~ 0.5, has 
been deduced for molecular clouds of a wide range of sizes (Larson 1981; Leung, Kutner 
& Mead 1982; Myers 1983; Dame et al. 1986; Myers k Goodman 1988; Falgarone, 
Puget & Perault 1992; Miesch &: Bally 1994). A turbulence spectrum characterized by 
m = 2 was also found in a H// region by Roy & Joncas (1985). It is interesting that 
this spectrum seems to prevail in different components of the ISM. The common factor 
in these cases is that the medium is compressible. Indeed, numerical simulations, with 
and without a magnetic field (Passot, Pouquet, & Woodward 1988, Vazquez-Samedani, 
Ballesteros-Paredes, & Rodriguez 1997), indicate that such a spectrum is characteristic 
to turbulence in a compressible medium (whether or not magnetic field fluctuations 
are important). 

The Kolmogorov spectrum strictly applies to the inertial wavenumber range for 
homogeneous and isotropic incompressible turbulence. In this range there is no energy 
input from the source generating the turbulence (input is at larger scales) nor energy 
losses (dissipation is at smaller scales). Thus, turbulent kinetic energy is cascaded from 
larger to smaller scales due to the nonlinear eddy interactions, at a constant rate that 
is independent of wavenumber. 

The present steeper spectrum, with m ~ 2, indicates that at each wavenumber a 
part of the cascaded turbulent kinetic energy is transformed to other energy forms. As 
a result, the turbulent kinetic energy rate is smaller the larger the wavenumber. In a 
compressible medium, part of the kinetic turbulent energy is indeed transformed to 
energy associated with density fluctuations. 

Finally, we note that the emitting neutral regions themselves, modeled here as 
thin slabs, are probably transient features created and disrupted by the large scale 
turbulence (Ballesteros-Paredes, Vazquez-Samedani, & Scalo 1999). Their life-times 
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are of the order of the timescales of the largest eddies. The spatial scales for which 
the power spectra were obtained are smaller than the above, and therefore the 
corresponding turbulence timescales are shorter than those of the very largest scales. 
With respect to the observed scales, the slab can be considered as a quasi-stationary 
structure in which the observed smaller scale eddies generate the density fluctuations 
that give rise to the intensity fluctuations. 

5.2. Turbulence Scales and Energy Sources 

The scale of the large-scale turbulence in the Galaxy should exceed 200 pc - 
the largest scale probed by the power spectra. The power-law shape of the spectra 
suggests that energy input to the turbulence is on yet larger scales. A probable energy 
source on these larger scales is the galactic differential rotation, as also suggested by 
numerical simulations of Wada & Norman (1999). In this natural scale will be 

a fraction of the galactocentric distance. However, for a 3- dimensional turbulence, the 
largest scale is likely to be comparable to the width of the galactic H77 layer of ~ 2kpc 
(Reynolds 1989; Reynolds 1991). 

In the SMC, the turbulence can be generated by a host of instabilities related to 
the bulk flows that result from the tidal interaction with the Galaxy and with the 
LMC, as well as to the expanding shells and super shells. The existence of large scale 
correlated density fluctuations in the SMC indicates the existence of turbulence on a 
scale ^ 4 kpc which is comparable to the size of the SMC . This observational evidence 
by itself is a very interesting result as the turbulence scale is an order of magnitude 
larger than that of the turbulence in the LMC, studied by Spicker & Feitzinger (1988) 
on the basis of radial velocities of 21-cm emission. 

The most plausible source feeding energy to such an extremely large scale turbulence 
are bulk flows induced by the tidal interactions with the Galaxy and the LMC, notably 
the very close encounter with the LMC that peaked about 2 x 10 8 yr ago and lasted 
for about 2 x 10 8 yr (Gardiner & Noguchi 1996). Such large scale bulk flows generate 
turbulence due to shear instabilities e. g. the Kelvin-Helmholtz instability (KH) or 
shock induced Richtmier- Meshkov instability (RM). These two mechanisms were shown 
to induce large scale turbulence in the intracluster gas of merging galaxy clusters by 
Roettiger et al. (1998) and by Goldman (1997), respectively. 

The simulations of Gardiner & Noguchi (1996) exhibit large scale (few kpc) velocity 
gradients with shear values S ~ 5 — lOkm/s/kpc, consistent with the (long period, 
population I) Cepheid radial velocity data of Mathewson, Ford,& Visvanathan (1986, 
1988). One expects that the large scale shear will produce turbulence with the largest 
eddies size It comparable to these spatial scales, and root mean squared turbulent 
velocities vt <; Sit = 10 4- 15km/s. This value is comparable to the sound speed of 
the embedding ionized medium, and supersonic with respect to the H/ regions, thus 
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enabling the density fluctuations in the latter. On smaller scales of ^ 1 kpc, also 
the giant super-shells observed in the SMC, contribute, to the turbulence. The giant 
supersehlls themselves are probably the result of a star formation burst that followed 
the interaction with the LMC (Zaritsky et al. 2000). 

The time scale for decay of the turbulence td ~ k/vt with a rms turbulent velocity 
of vt ~ lOkm/s and scale k = lkpc is ~ 4 x 10 8 yr. Thus, we are witnessing a turbulence 
that didn't have enough time to decay, since the close encounter with the LMC of 
~ 2 x 10 8 yr ago. Moreover, the simulations of Roettiger et al. (1998) and simulations 
of RM instability (Rotman, 1991) indicate that continued excitation (even at a lower 
level than the original one) can keep such a turbulence going for times considerably 
longer than the above estimate. Such a lower level excitation can be provided by 
ongoing tidal interactions of smaller' strength. 

Another question of interest is the near constancy of the slope of the spectra. We 
know that small scale (from lkpc and down) ) energy sources exist (supergiant shells, 
smaller scale supernovae remnants, smaller scale hydrodynamic instabilities) as well 
as smaller scale energy sinks (e.g. magnetic excitations). Consider a source /sink a 
at wavenumber k in an interval A A; such that k ^ Ak >> ko with ko corresponding 
to the largest scale in the turbulence. Denote by e and Ae the energy rate due to the 
large scale turbulence and the change due to the small scale source, respectively. Since 
F(k) oc e(fc) 2/3 the jump in the turbulence spectral function 

„ 2 Ae 

F(k) 3 e { ' 

Even for a relative change of ~ 20%, the jump will be <J 10%, consistent with the 
scatter of the observational spectra. 



5.3. Implications of the Turbulence on the SMC Magnetic Field and 

Abundances Gradients 

The consequences of this turbulence may be quite important regarding magnification 
of large scale magnetic fields, and enhancing diffusion coefficients. In the SMC, the 
turbulence largest eddies are comparable in scale to the size of the SMC. This implies 
that turbulent mixing should smooth any large scale abundance gradients. This 
seems indeed to be the case observationally (Dufour 1975; Kobulnicky 1998) . The 
turbulence is also expected to amplify and shape up the large scale magnetic field. 
Indeed, the SMC observational data indicate the existence of a large scale disordered 
field (Haynes et al. 1991; Ye & Turtle 1991). The observed magnetic field is ~ 6.5/iG. 
a value close to that expected from a dynamo amplification to an equilibrium value 
B = 5/iG(ft/10kms _1 )(n/lcm~ 3 ) 1 / 2 with vt the turbulent root mean squared velocity 
and n the number density in the ionized embedding region. 
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5.4. Concluding Remarks 

The interpretation of the Hj power spectrum of the SMC in terms of an underlying 
large scale turbulence seems natural and as noted above is consistent with other 
observational data. The existence of a large scale turbulence in the SMC gains also 
support from the recent finding of Stanimirovic et al. (2000), that the autocorrelation of 
the dust column density is also a power law with an index ~ —3. In our interpretation, 
the dust will be just another passive marker waved by the turbulence and thus should 
have the same power-law spectral index as the Hj intensity. 

Moreover, observations of stars of various ages were used by Zaritsky et al. (2000) 
and by Kunkel, Demers, & Irwinl (2000), to conclude that hydrodynamic interactions, 
following the recent close passage of the LMC, rather than just gravity are responsible 
for the morphology and dynamics of the ISM of the SMC. 

Comments by the referee and by Dr. Steven Shore have contributed to the 
improvement of the presentation. This work has been supported by the US-Israel BSF 
grant 94-314 and Israel Science Foundation grant 561-9911.2. 
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